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Abstract 

Ultrahigh energy neutrinos can provide important information about the distant 
astronomical objects and the origin of the Universe. Precise knowledge about their 
interactions and production rates is essential for estimating background, expected 
fluxes and detection probabilities. In this paper we review the applications of the high 
energy QCD to the calculations of the interaction cross sections of the neutrinos. We 
also study the production of the ultrahigh energy neutrinos in the atmosphere due to 
the charm and beauty decays. 



1 Introduction 



The ultrahigh energy neutrino physics has attracted a lot of attention during last couple of 
years. Neutrinos are important messengers of information on distant extragalactic sources 
in the Universe and the neutrino astronomy has many advantages over the conventional 
photon astronomy. First of all, neutrinos, unlike photons, interact only weakly, so they 
can travel long distances being practically undisturbed. The typical interaction lengths 
for the neutrino and the photon at energy E ~ 1 TeV are about 

L\ nt ~ 250 x 10 9 g/cm 2 , CJ nt ~ 100 g/cm 2 . 

Thus very energetic photons with energy bigger than 10 TeV cannot reach the Earth from 
the very distant corners of our Universe without being rescattered. On the other hand 
neutrinos can travel very long distances and are also not deflected by galactic magnetic 
fields. 

At ultrahigh energies the angular distortion of the neutrino is very small. Typically 
the angle between the neutrino and the produced muon is about 

H ~ (£„/TeV)°-7 " 

Therefore the highly energetic neutrinos point back to their sources. The interest in the 
neutrinos at these high energies lead to the development of the experimental devices such 
as AMANDA 1 , NESTOR [2j, ANTARES or planned ICE CUBE gj observatories. 
For the reliable observation of the neutrinos, their cross sections and production rates 
in hadronic matter have to be well known. Even though the neutrinos interact only 
weakly with other particles, strong interactions play an essential role in the calculations 
of the neutrino production rates and their interaction cross section. This is due to the 
fact that neutrinos are coming from the decays of various mesons such as tt,K,D and 
even B which are produced in the high energy proton-proton (or proton- nucleus, nucleus- 
nucleus) collisions. These hadronic processes occur mainly in the atmosphere though 
possibly can be also present in the accretion disc in the remote Active Galactic Nuclei. 
Also, the interactions of highly energetic neutrinos with matter are dominated by the 
deep inelastic cross section with the nucleons/nuclei as described in Sec. 4. This is why 
the knowledge about QCD from high energy collider experiments such as HERA and 
TEVATRON is invaluable. The production and interaction of highly energetic neutrinos 
with nucleons involves gluon distribution probed at very small values of Bjorken variable 
x (~ 10~ 9 ) which corresponds to a very high c.m.s. energy. This is the kinematic domain 
which is currently inaccessible in these colliders (at HERA the minimal x m i n ~ 10~ 5 at 
Q 2 ~ lGeV 2 , where — Q 2 is the virtuality of the photonic probe), and this calls for a 
reliable extrapolation of parton densities constrained by the current experimental data. 
In this review we concentrate mainly on the application of the knowledge and experience 
gained in testing high energy QCD at present colliders to the problem of precise evaluation 
of the neutrino-nucleon (or nucleus) scattering cross sections. In addition, the production 
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rate for the atmospheric neutrino flux from the charm or beauty decays (the so called 
prompt neutrino flux) will be analyzed. 

This paper is organised as follows: in the next section we briefly discuss the origin of the 
ultrahigh energetic neutrinos. In Sec. 3 we describe the mechanism of the prompt neutrino 
production at high energies in the atmosphere and present the calculation of the cross 
section using different models for the small x dynamics, including parton shadowing. We 
then follow the sequence of fragmentation, interaction and decays of the charmed mesons 
to obtain the neutrino fluxes at high energies. In Sec. 4 we discuss the interaction of 
the neutrinos with the nucleons and in particular we evaluate the deep inelastic neutrino- 
nucleon cross section based on the various parametrisations of the parton distribution 
functions. In Sec. 5 we follow the process of the neutrino transport through the Earth 
and present the angular dependence of the fluxes after travelling through the matter. Also 
in Sec. 5.1 we discuss the effect of the v T regeneration. Finally, in Sec. 6 we discuss the 
effects of the neutrino oscillations. 



2 Sources of high energy neutrinos 

There are various sources of ultrahigh energy neutrino fluxes. The currently measured 
ultrahigh energy neutrino flux is the atmospheric one which results from interactions 
of cosmic rays with the nuclei in the atmosphere. It is rather steeply falling flux E~ 3 ' 7 , 
however it will have the large energy tail due to the charm production which will be 
discussed in detail in the next section. 

There are several predictions for the ultrahigh energy neutrino fluxes other than the 
atmospheric one, most of them differ by orders of magnitude. For the detailed review on 
the possible ultrahigh neutrino sources see for example E] arid the recent constraints 
on the ultrahigh energy neutrino fluxes have been given in [7j. 

Here we briefly review the possible sources of ultrahigh energy neutrinos. 

One can expect a contribution to the galactic flux resulting from interactions of cosmic 
rays with the interstellar medium, and the extragalactic flux which originates from the 
pp interactions of cosmic rays with intracluster gas in the galaxy clusters. At the very 
high energies one can also expect cosmogenic neutrino flux coming from the interactions 
of the highest energy cosmic rays with the ubiquitous photons from cosmic microwave 
background jH]. 

An interesting possible source of the ultrahigh energy neutrinos are the Active Galactic 
Nuclei |Hl llUllTTlll2j . These are highly luminous objects, which luminosity from the nucleus 
reaches or exceeds that of our Milky Way galaxy Lagn — 0.1 — lOLMiikyWay The energy 
of the AGN comes from gravitational energy of a supermassive black hole in the center 
of the nucleus with mass Mbh = 10 6 — 10 10 Af sun . The matter accrets onto a disc which 
surrounds a black hole, and the hot gas from inner parts of the disc can move away from 
it and get squeezed by the magnetic field. This leads to the emergence of highly luminous 
jets which can extend over the distance of 1 — 50kpc. There are different classes of AGNs 
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like Seyfert, blasars and quasars. Most probably they are similar objects but with just a 
different orientation towards the Earth. The production of the neutrinos occurs through 
the reaction of pj — > A + — > 7r + + n with subsequent decay of the pion. Protons can 
be accelerated in the accretion disc and interact there with the photons, or in the jets 
where they can interact with the synchrotron radiation photons. The shock acceleration 
has been proposed as a most plausible mechanism for the acceleration in AGNs (see for 
example 

Gamma Ray Bursts which are characterised by the intense and short-time photon 
spectra can also be considered as a possible source of neutrino flux. The detailed origin 
of GRB is not known, though most probably they are very massive hypernovae |13| . The 
mechanism of GRB is modelled via expanding relativistic fireball. If the protons are also 
shock accelerated they might interact with the synchrotron photons to produce neutrinos 
via pion photo-production |14j . 

Super novae are powerful sources of the neutrinos and gamma-rays at nuclear energies. 
Though, there has been no evidence for the high-energy gamma ray or neutrinos fluxes 
from supernova SN1987A, it is possible that supernovae could be sources of the highly 
energetic neutrinos. 

Other, more exotic mechanisms are also studied. In particular the annihilation of 
Weakly Interacting Particles (WIMPs) trapped in the cores of Sun or the Earth |15j . 
WIMPs have been discussed in the context of the dark matter component of the universe. 
Using the predicted density for the dark matter one can estimate the flux of neutrinos 
which would come from the process of WIMPs annihilation. For example within the 
supersymmetric scenario annihilation of neutralinos has been proposed as a source for 
the flux of muon neutrinos |15j . In the supersymmetric extension of the standard model 
with conserved R-parity neutralino is the lightest SUSY particle-LSP and it is stable. 
Neutralinos could be therefore attractive candidates for the dark matter WIMPs. If their 
density is high enough they could annihilate producing ultrahigh energy neutrino flux in 
the final state ^5] (for recent calculation see |16j^. 

Another speculative source are so called top-down models j^l El which assume the 
existence of the very heavy dark matter particles I2()j - Wimpzillas. In various GUT 
models these heavy particles occur naturally with a mass around the GUT scale mx — 
2 x 10 16 GeV. If their lifetime is long enough they would contribute to the cosmic and 
neutrino flux at very high energies. In some models very heavy particles are produced 
continuously from the emission of the topological defects |21l 1221 123j that are relics of the 
cosmological phase transitions in the early universe. For example, during the reheating 
phase in the early universe the ultraheavy relics can be produced quite copiously and they 
could be the candidates for the dark matter in the universe. These supermassive particles 
with the relatively long lifetime could be concentrated in the galactic halo |24| . This is 
very interesting possibility since it could explain the possible excess of the data [2E] beyond 
the GZK j2H] cutoff: the products of the decays of these wimpzillas would have energies 
exceeding 10 19 eV. This in consequence would also lead to the enhancement neutrino flux 
at the ultrahigh energies. 
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Finally, we stress that an important general upper bound on the ultrahigh energy 
neutrino flux has been derived [27j from the known cosmic ray flux. This limit is applicable 
to the AGNs and GRBs and could only be exceeded through the exotic scenario (like 
WIMPs annihilation or the top-down model ) or if the production source for the cosmic 
rays is optically thick to proton-photon pion production or proton-nucleon interactions (so 
called "neutrino only factories"). 



3 Atmospheric flux and prompt neutrino production 

In this section we study in detail the production of the high-energy neutrinos in the 
atmosphere. This neutrino flux constitutes an important background to the flux of cosmic 
neutrinos discussed in the previous section. The atmospheric neutrino flux is produced by 
the interaction of cosmic rays with the nuclei in the atmosphere |28| . The neutrinos mainly 
come from the decays of pions and kaons. At energies of around 100 TeV the decay length 
of these mesons is so large, that they interact before decaying and loose substantial amount 
of energy. This effect leads to the steeply falling spectrum of atmospheric neutrinos [251 
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At higher energies E v > 100 TeV the atmospheric neutrino flux will be dominated by the 
prompt neutrino production originating from charm and also beauty hadron decays. These 
hadrons (D , D+(~\ D+ mesons, barion A c and their B counterparts) are characterised 
by a very short lifetime ~ 10~ 12 s, so up to energies 10 9 GeV they will not interact before 
decaying and the flat flux of the neutrinos will be produced which can extend up to these 
energies. This prompt neutrino flux at high energies is a background for the sought after 
cosmic neutrinos and its precise evaluation is essential. 

3.1 cc production cross section 

There are several QCD-motivated calculations of the prompt neutrino flux |30l 1311 132] . 
Below, we present a particular calculation which includes small x resummation effects as 
well as estimate of the gluon shadowing [3*3] . 

The cross section for the charm production can be evaluated within perturbative QCD 
framework. The high energy regime of the production process pp — > cc + X probed by the 
cosmic rays is however not accessible in the present colliders and reliable extrapolations 
of these predictions are therefore necessary. 

In Fig. n we present the diagram for the charm production in the collision of two 
protons, x F = E c /E is the Feynman variable characterising the fraction of the energy of 
the incident cosmic ray proton carried by the produced charm quark. The cross section 
for the production of the cc pair in the proton-proton collision reads 

da pp ^ c+x f do -^ 5 c 

dx x dx 2 dz g(xx,fi F ) — £ f zl££ g{x 2 , fi F ) S(zxi - x F ), (1) 



dx F J dz 
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Figure 1: QCD diagram for the prompt charm production in the proton-proton collisions. 



where ° 9 ^ cc is the partonic cross section for the process gg — ► cc and g(x,fi 2 ) is the 
gluon density in the proton evaluated at the momentum fraction x and the scale /i 2 . In 
the domain of high energies that we are interested, the typical momentum fractions are of 
the order x± ~ xf ~ 0.2 — 0.1 and consequently X2 — M 2 5 /2xfs ~ 10 -4 — 10~ 9 where M 2 5 
is the mass of the cc system. In general the higher the energy the lower the momentum 
fraction X2 in question. It means that in this process the gluon density g(x2, /J- 2 ) in Eq. Q 
is probed in the regime of X2 which is outside the range of any present accelerator data. 
We consider three approaches for an extrapolation into the low x region: 



• extrapolation of the standard DGLAP evolution 

• united BFKL/DGLAP evolution equations 

• inclusion of saturation effects, Golec-Biernat and Wuesthoff model 



We briefly describe these calculations below. The most straightforward approach is 
to use the standard gluon distribution obtained within the DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) framework from the global fit [HI] and extend it to lower values 
of x using the double leading-log approximation 



g(x,Q 2 ) = ,(x = 10-Q 2 )exp Jl^ln^-ln^- 

V V b as{Qo) x 



16Nc ln as(Q) ln 10- 5 



as(Qo) x 

where nj = 4 and b = 25/3. The running leading-order (LO) coupling is defined as 

47T 



cts{Q) 



&log(QVA| CD ) 
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This approach correctly incorporates the leading a s lnQ 2 lnl/x terms. In this calculation 
(called MRST on the plots) we have used the LO formulae for the partonic cross section 
and multiplied them by the factor K = 2.3 which accounts for the next-to-leading (NLO) 
corrections. 

Another approach is to use the BFKL(Balitsky-Fadin-Kuraev-Lipatov) |35| framework 
which resums the leading terms a s \ul/x. This approach is probably more appropriate 
since in the process considered, Q 2 ~ m 2 is fixed and not too large whereas the values of 
Bjorken x probed are very small. We have chosen to use the gluon distribution function 
obtained from the solution of unifed BFKL/DGLAP evolution equations |36| . This ap- 
proach (which we call KMS on the plots) has the advantage of treating both the BFKL 
and DGLAP evolution schemes on equal footing and it also resums the major part of 
the NLO corrections to the BFKL equation via the so-called kinematical constraint |37j . 
In this framework one is using the high energy factorisation theorem |38j together with 
the unintegrated gluon distribution function f{x,k 2 ) which is related to the integrated 
distribution xg(x,[i 2 ) 

rf? dk 2 

xg(x,n 2 ) = / ~^rf(x,k 2 ) . 

Here k 2 is the virtuality of the gluon. More precisely for the evaluation of the cross section 
one has to use the off-shell matrix element for gg — > cc convoluted with functions f(x, k 2 ) 

da PP ~* C+X f ^.off-shell 

dx\ dx2 dz dk\ dk\ f(x\, k\) f{x2, k 2 )5(zxi — xf) , (2) 



dxF J dz 

where kf and k\ are the virtualities of the two incoming gluons and da ° S d ^ cU is an off-shell 
matrix element for the gg — ► cc process. This approach is more appropriate at small x 
since it accounts for the resummation of the a s In 1 /x terms in the splitting gg — > cc. 

In the third calculation (which we refer as to GBW) we have included the effects which 
might come from the parton saturation. In the region of high energies, where the parton 
densities become high, the gluon recombination effects can become important. These 
effects lead to the slower increase of the parton density with energy and in consequence 
to the damping of the cross section. This is called perturbative parton saturation. For a 
review on various approaches to this phenomenon, see [SHI- In our calculation we have used 
very successful model by Golec-Biernat and Wuesthoff 40 j which is formulated within the 
dipole picture of high energy scattering. This model has been successfully applied to the 
description of the inclusive and diffractive data on deep inelastic scattering at HERA. Let 
us briefly sketch the philosophy of this approach which was originally formulated for the 
deep inelastic process of the scattering of electron on a nucleon. In the dipole picture the 
virtual photon fluctuates into a qq dipole pair and then interacts with the target. These 
two processes are well separated in time. The total cross section for the "f*p scattering 
can be expressed in the following form 

tTT,L(x,Q 2 ) = E / d?rdz\& TL (r,z,Q 2 )\ 2 aW(r,x) , (3) 
q J 
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Figure 2: The Z moment of charm production in proton- nucleus collisions plotted as a 
function of energy E c of the produced charm quark. Solid line - GBW saturation model, 
dashed line - standard DGLAP calculation based on MRST partons, dotted line - KMS 
prediction based on the unifed BFKL/DGLAP evolution equations. The Z moment is 
shown for illustration. In the calculation of the neutrino fluxes the differential cross section 
is convoluted with the primary cosmic ray flux. 



where Vl/Ip L is the virtual photon wave function (T stands for transversely and L for 
longitudinally polarized virtual photons) and a q ^(r, x) is the qq dipole cross section which 
contains the information about the interaction between the qq dipole and the proton. 

The r is the dipole pair transverse size, and z is the longitudinal momentum fraction 
of the photon carried by the quark. The photon wave functions $m £ are well known 
in leading order. In the GBW saturation model 0U] the dipole cross section has been 
parametrised in the following form 



1 — exp(- 



r Q s (x). 



(4) 



with the saturation scale 



Qs(x) 



Qo[- 



-A/2 



(5) 
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which governs the x (or energy) behaviour. The parameters of the GBW model are 
Qo = 1 GeV, xo = 3 x 10~ 4 , A = 0.28, uq = 23 mb and are chosen to give the best 
description of the experimental data at HERA collider. 

This model has the property that for small dipole sizes r < 1/Q s (i) the cross section 
is also small and behaves proportional to r 2 . For large dipole sizes r S> 1/Q s {x) the dipole 
cross section a qq does not grow anymore, but rather saturates to the constant value given 
by parameter uq. The transition scale Q s (%) between the linear and saturated regimes 
called saturation scale grows with energy like a power (|5|). The larger the energy the 
smaller the dipole which will saturate. This model has been extended from 7*5 — > qq to 
describe the gg — ► cc process as well and applied to the prompt charm production [SHj. Let 
us finally note, that GBW model contains much more than the perturbative saturation. 
This model is valid for all values of Q 2 and the dipole cross section extends to arbitrary 
large dipole sizes, therefore GBW parametrisation necessarily embodies a lot of soft or 
nonperturbative components. 

In Fig. 12 we have summarised the results for the charm production cross section based 
on the three above-mentioned calculations. We have plotted a moment 



he calculation labelled MRST is the standard DGLAP approximation, KMS includes ad- 
ditional small x dynamics through the unifed BFKL/DGLAP evolution whereas GBW ac- 
counts for the gluon saturation. We see that, the calculation based on unifed DGLAP/BFKL 
equation gives the highest prediction for the energies > 10 6 GeV, whereas the calculation 
which contains the gluon saturation tends to have slowest increase with energy. This 
is because the parton distributions obtained from solution to the unifed DGLAP/BFKL 
evolution equations exhibit a stronger growth with energy than the ones from standard 
calculation within the DGLAP approach. The difference between the KMS and MRST 
calculations at lowest energies can be explained by the fact that the latter has been multi- 
plied uniformly by the next-to- leading K factor (~ 2.3), whereas the KMS calculation uses 
the off-shell matrix element which in principle contains subleading corrections, important 
at high energies. On the other hand, in the low energy regime the gluon distribution is 
probed at rather large values of x ~ 0.1 and so the In 1/x resummation in matrix element 
in KMS calculation is not effective. The KMS prediction should be in principle multiplied 
by the energy dependent K factor which would vanish in the high energy region. 

3.2 Calculation of the prompt and v T flux. 

In order to calculate the final neutrino flux one has to follow the sequence of the charm 
production, charm quark fragmentation and charm meson decay into neutrinos. Addition- 
ally the attenuation of the cosmic rays has to be taken into account as well as possibility 
of the interaction of the charmed mesons. This sequence of processes has been followed 
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Figure 3: The prompt v^ + v^ and v T + v T fluxes, arising from cc production, fragmentation 
and decay, obtained using the three different extrapolations of the gluon to very small x. 

by use of transport equations (see for example ) which are differential equations in 
depth X of atmosphere traversed by the particle 

X = p{ti) dti , 
Jh 

where p(h) is the density of the atmosphere as measured at height h. The density profile 
has been taken to be exponential p(h) = po exp(— h/ho) with po = 1.6 x 10~ 2 g/cm 3 and 
ho = 6.4 km. The transport equations which determine the fluxes have the following form 
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Figure 4: The prompt (a) + and (b) v T + P T fluxes calculated using the charm 
production cross sections obtained from GBW model. Also shown are the conventional 
muon and tau neutrino atmospheric fluxes, where the latter originates, via neutrino mixing 
transitions, from the former. There is also a contribution to the prompt v T + v r flux from 
beauty production, which is not included here, but is shown in Fig. 5. 



d<f> N (E,X) 



1 



dX An(E) 

Me,x) = 

dME,X) 1 



n(E, X) 

oo 



dE'dx c ct> N {E',X) — ^- 5(E-x c E'), 
E P dx c 



dX 
dHE,X) 
dX 



f°° dn 



/ dE'dx<t)i{E',X) 



1 



Xf QC (E',X) 



dx 



(6) 



Here, 4>n(E,X) is the cosmic ray flux at the depth height X, with 4>n(E,0) being the 
initial flux (primary cosmic ray flux). (j)c(E, X) is the c quark flux, (f>i(E, X) is the charmed 
hadron flux (where i = D°, D^, D°, D^,A C ) and <fti(E, X) is the lepton flux. The nucleon 
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Figure 5: The prompt v T + v T fluxes originating from cc and bb production and decay, 
which respectively arise from the D s — ► tv t decay, and from the B^, B°, B s and A& 
semileptonic tv t decay modes. The upper plot shows the breakdown into the direct v T 
contribution (continuous curves) and the indirect r — > u T contribution (dashed curves). 
The lower plot shows the total prompt v T + U T flux, together with its components of cc 
and bb origin. Also shown is the non-prompt v T + v T flux arising from — * v T oscillations 
from the conventional atmospheric flux. 



attenuation length is defined as 



An(E) 



1 — Zmn 



where Znn is the spectrum- weighted moment for the nucleon regeneration and A at is the 
interaction thickness. 

The attenuation length for charmed hadrons has two terms: decay length and interac- 
tion length 

1 1 1 - Z c - C 



+ 



A; Af c Xi ' 

where A^ ec is the decay length and Aj/(1 — Z cS ) is the attenuation length for the charmed 
hadron with Z c5 being the charm regeneration factor. The branching fraction of the decay 
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of the charmed hadron into leptons is given by B{i — > I). The differential cross section 
a J x ~* c is given by the calculation presented in the previous section. The distribution dn ^ x 1 
describes the fragmentation of the charm quark into the hadron and we have taken it to be 
proportional to 5(x£> — 0.75x c ) rather than using complete fragmentation function. Finally 
the lepton distribution from the hadron decay n dx has been taken from PYTHIA Monte 
Carlo HI]. 

In Fig. El we present the predictions for the and v T fluxes coming from the solution 
to the transport equations © based on three models of the gluon extrapolation in — 
to small values of x. The KMS prediction is highest at high energies due to the steepest 
increase of the gluon density caused by the resummation of the a s \n{l/x) terms. On 
the other hand the GBW saturation model gives lowest predictions, due to the gluon 
saturation effects embodied in the model. To be more precise we have considered several 
scenarios for the saturation(for details see [SB]); see Fig. 21 where the band in both plots 
shows the uncertainty of the model due to the different assumptions (like the strength of 
the 3-pomeron coupling). The prediction for the Vn flux is compared with the standard 
atmospheric flux coming from the decay of lighter mesons: it and K. The prompt v T 
flux is practically the only source for the atmospheric neutrinos at energies E > 10 TeV. 
Additionally we have a v T component which comes from the oscillations — > u T , dominant 
at lower energies E < 1 TeV. 

Apart from charm decay it is also possible to have neutrino fluxes from the beauty 
hadron decays. For the f^,f e fluxes these processes are negligible, but for the v T the B 
decays are an important source of the flux. This happens because due to larger mass of 
B hadrons with respect to charmed hadrons more decay channels to v T are opened (in 
the case of charm only D s gives contribution to u T ) and this compensates for the lower 
value of the production cross section. In Fig. \5\ we show the predictions for the v T flux 
from bb decays as compared to cc decays and we see that at energies around 10 6 GeV 
the bb component constitutes about 30% of the total v T flux, and at higher energies this 
proportion is even bigger. 



4 Neutrino cross sections at high energies 

Another topic that relates QCD at small x and the ultrahigh energy neutrino physics is 
the problem of the neutrino interactions with the encountered matter on their way to the 
detector. The highly energetic neutrinos will interact with electrons and nucleons in the 
Earth and the detailed knowledge about the cross sections is essential for the knowledge of 
the rate of muons to be observed in the detector. There are various types of interactions 
of neutrinos. The interactions with electrons of the muon and electron neutrinos can be 
listed as follows 

and the interactions of the electron antineutrinos 

v e e — ► D e e, D e e — > z^/i, v e e — > hadrons . 
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Figure 6: The total vN charged current cross section and its decomposition into com- 
ponents of different origin (valence, sea quarks) as a function of the laboratory neutrino 
energy E v . 



The first group of interactions is characterised by the constant cross section at high ener- 
gies. On the other hand the interactions u e e — > anything have a characteristic sharp peak 
coming from a Glashow resonance at E u ~ 6 PeV due to the W~ production. 

This window of the enhanced electron antineutrino interactions together with the so- 
called double bang events in the case of tau neutrinos [32j can be used to disentangle the 
flavour composition of the incoming cosmic neutrino flux in the neutrino observatories. 

The interactions with nucleons proceed through the charged and neutral current deep 

CC NC 

inelastic process vN — ► IN, uN — ► vN. The cross section for the charged current in- 
teraction, vN — > IN scattering on the isoscalar target iV = has the following form 



d 2 a cc 2G\M N E V ( M? 



dxdy n \Q 2 + M? 



w 



w 



\ ■ [xq(x,Q 2 )+xq(q,Q 2 )(l-y) 2 ] , (7) 
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Figure 7: The total vN cross section and its decomposition into charged current and 
neutral current components as a function of the laboratory neutrino energy E v . 

where Mn mass of the target nucleon, Q 2 is the exchanged vector boson four momentum 
transfer, x = Q 2 /(2Mu) Bjorken variable, y = v/E inelasticity (where v = E v — E\ is 
energy loss in the lab frame), Gp Fermi coupling constant and My/ mass of the W boson. 
xq(x,Q 2 ) and xq{x,Q 2 ) are the quark (antiquark) momentum distributions evaluated at 
the momentum fraction x and at scale Q 2 . 

The characteristic feature of this deep inelastic cross section is the fact that the parton 
densities grow with decreasing values of x as xq(x,Q 2 ) ~ x~ x as we know from the 
measurements of the proton structure function at HERA collider [13]. This results in the 
strong growth of the total vN cross section with energy 

a tot {E) = / dxdy^- ~ E x . (8) 



dxdy 

The behaviour of the charged current cross section with energy is illustrated in Fig.Elwhere 
the QCD calculation 44 of the partonic densities has been done within the framework of 
the unifed BFKL/DGLAP equations [HBj. 
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Figure 8: The contour plot showing the (x,Q 2 ) domain of the dominant contribution 
to the din(i/x)diii(Q' 2 /A' 2 ) f° r total vN interaction for neutrino laboratory energy E = 
10 11 GeV. The 20 contours are such that they enclose contributions of 5, 10, 15, ... % of 
the above differential cross section. 

We show various components to the cross section. First thing to note is the fact that 
above 10 5 GeV sea quarks(ucf, cs) start to dominate over the valence distribution. The 
valence quark contribution grows linearly with energy for E < 10 5 GeV due to the opening 
of the phase space, and then flattens out because the factor M^/(Q 2 + My^) provides a 
cutoff and because of the fact that the number of the valence quarks is finite. On the other 
hand the density of sea quarks grows as a power with x which is translated into a power 
growth with energy as in (JHJ) . Threshold effects in the production of heavy quarks ( bt ) 
are also visible at lower energies. In Fig. Q we show the total cross section together with 
the charged and neutral current components. The neutral current cross section constitutes 
about 30% of the total cross section at high energy, the fraction being determined by the 
ratio of combinations of the chiral couplings. 

One of the important questions is to what extent are these predictions stable. After 



16 



all, one is extrapolating these calculations a long way into the new kinematical domain. 
There are various sets of the parton distributions which can have different asymptotic 
behaviour at very small x. However, as has been investigated in |44| and |45j the parton 
distributions are quite well constrained by the current experimental data and so variations 
of these distributions result at most in factor 2 — 3 difference in normalisation at highest 
energies. 

To visualise the kinematic regime probed in ultrahigh energy neutrino-nucleon in- 
teractions we have made the contour plot, see Fig. |HJ of the differential cross section 
dini/xd\nQ i /A i m ( x > Q 2 ) plane. The contours enclose the regions with different contribu- 
tions to the total cross section a{E u ). We see that for the very high energy E = 10 11 GeV 
the dominant contribution comes from the domain Q 2 ~ Mj^ and x m \ a ~ M^/(2MnE) ~ 
1CP 7 where Mjv is the nucleon mass. 



4.1 High partonic densities and saturation 

In the limit of very high energies one might expect the shadowing (or saturation) correc- 
tions to appear and modify the asymptotic behaviour of the total cross section. In Ref. |33| 
it was argued, on the basis of the GBW saturation model [IP] , that due to the fact that the 
saturation scale is typically much smaller than the average value of the momentum transfer 
Q s {x) <C {Q) ^ M\y, the actual size of the shadowing corrections is negligible, since the 
contribution to the cross section from the nonlinear regime Q < Q s (x) is very small. A 
more detailed estimate has been performed 07j where the nonlinear Balitsky-Kovchegov 
|48j has been solved which takes into account the saturation effects. It was shown that, 
even though the typical values of momentum transfer Q 2 are in the linear regime, the 
nonlinear effects can tame the growth of the total cross section and reduce the magnitude 
at highest energies ~ 10 12 GeV by a factor of 2 or so. It was shown in 023J that the 
partonic distribution persists to feel the nonlinearity even for momenta higher than the 
saturation scale, more precisely in the window Q^(x)/ Aq CD > Q 2 > Q 2 {x). Therefore 
the saturation effects can become non-negligible even at highest energies and they could 
be enhanced by the scattering of the neutrino on the nucleus. 



5 Penetration through the Earth 



We have mentioned at the beginning that the neutrinos interact very weakly and therefore 
do not get disturbed during their way from the source to the detector on Earth. However, 
for the ultrahigh energy neutrinos the situation changes due to the increasing value of 
the neutrino-nucleon cross section, see Figs. 16171 The ultrahigh energy neutrinos will 
undergo attenuation when penetrating Earth j45j . It is well illustrated in Fig. |§] where the 
interaction length L^{E y ) = ct „jv^ ^ Na is shown together with the maximal column depth 

^Earth = 1-1 x 10 10 cmwe which is encountered by neutrino traversing along the diameter 
of the Earth. It is obvious that the interaction length becomes shorter for neutrinos with 
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9) 




E v [GeV] 



Figure 9: The interaction length for the neutrino interaction with isoscalar nucleon target 
obtained from the KMS model. The straight line corresponds to the maximal column 
depth length which can be encountered by the neutrino traversing the Earth through the 
middle. 

energies exceeding 40TeV. Thus Earth will start to become opaque to highly energetic 
neutrinos 



These effects can be studied using transport equations for the neutrino flux 
dI(E,T) , f dy da NC {E',y) 



dr 



-v TOT (E)I(E,t) + / *V- *>»C{E ,y) 
J 1 - V ay 



where the total cross section is the sum of charged and neutral current components 
otot = cr^vc + a cc an d where y is the fractional energy loss E' = r = NaX(z) = 
Na Jq dz'p(z') with p(z) being the density of the Earth along the neutrino path and Na 
is Avogadro number. 

The total depth length r depends on the incident nadir angle which is defined as the 
angle between the direction of the neutrino arrival and the normal to the Earth surface 
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passing through the detector. The first term in Eq.© corresponds to the neutrino atten- 
uation whereas the second term is responsible for the neutrino regeneration through the 
neutral current interactions at higher energies. It is useful to study the general properties 
of the solution by showing the shadowing factor S(E,r) = I(E,t)/Iq(E) for two different 
forms of initial flux, Fig. I1UI First is the steeply falling atmospheric |29j neutrino flux 
Io ~ oE -3 ' 7 and second is the Active Galactic Nuclei flux which is characterised by 
rather flat behaviour in the interval 10 3 < E < 10 5 GeV. 

Three different incident angles are considered in Fig. ^] which correspond to different 
values of X = 1.1 x 10 10 , 0.45 x 10 10 , 0.072 x 10 10 cmwe for 6 = 0°, 40°, 80° correspondingly. 
Additionally we show the pure attenuation factor A = exp(— <7tot"7") which would occur, if 
the neutral current regeneration given by the second term in Eq.@ would not be present. 
It is apparent that for the large column depths, the attenuation at high energies is very 
strong, essentially only about 10% of the initial flux survives at energies 100 TeV and the 
flux gets completely absorbed at yet higher energies. What is also visible is the fact that 
the regeneration due to the neutral current interaction is most prominent for flat fluxes like 
the one from AGN, and rather negligible for steeply falling spectra like the atmospheric 
one. In some range of energies the effect of regeneration is quite large, in which case the 
shadowing factor is bigger than 1, see Fig. 1101 

In Fig. ^2 we show the resulting neutrino flux 1(E) reaching the detector at different 
nadir angles 6 as a function of energy for GRB model [13], Top-Down model ^H] and the 
atmospheric flux. The GRB and TD fluxes dominate over the atmospheric flux at energies 
above 10 5 GeV. However at small nadir angles we see the strong attenuation of all fluxes 
at higher energies. The smaller the nadir angle the greater the shadowing. Therefore one 
has to adopt different method for searching the ultrahigh energy neutrino events in the 
detectors whenever E > 1 PeV. For smaller energies it is sufficient to look for the upgoing 
muon events, which are always neutrino induced. This reduces the background from the 
atmospheric muons. For higher energies due to the discussed Earth shadowing one has 
too look for the horizontal or downgoing highly energetic muon events. A detailed study 
on the possibility of detecting tau neutrinos using horizontal arrays has been performed 
in |5I]. 

5.1 Regeneration of r neutrinos 

It has been pointed in [52] and further studied in [53 [53 that the propagation of tau 
neutrinos through the Earth is quite different from the muon and electron neutrinos. The 
Earth in principle never becomes opaque to the high energy tau neutrinos [52] ■ This is 
due to the very short lifetime of the r lepton, which decays quickly, producing in the final 
state tau neutrino still with very high energy. This process continues until the v T energy is 
reduced to the value at which the interaction length is equal to the column depth through 
the Earth. The tau neutrinos will thus never disappear from the high energy spectrum, the 
flux will be only redistributed towards the lower energies with possible pile-up around the 
transparency energy. This is different than in the case of muons, which live substantially 
longer and loose their energy via interactions. When they finally decay, they produce a 
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Figure 10: The shadowing factor S(E, r) = -prngy as a function of the neutrino energy E v 
for three different incident angles 9 = 0°,40°,80° and two different models for the fluxes: 
atmospheric (left plots) and AGNs (right plots). The solid line corresponds to the solution 
of the full transport equation whereas the dashed line corresponds to the shadowing factor 
being purely absorption S(E,t) = ex.-p(—a tot (E)T) without regeneration. 
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Figure 11: The neutrino flux I(E,t) as a function of the neutrino energy E v for four 
different incident angles 9 = 0°,40°,80° and 180° (initial flux I°(E)) and three different 
models for the fluxes: atmospheric, TD model and the GRB model. 
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muon neutrino with very low energy. The effect of tau regeneration is mostly prominent 
for the flat fluxes, where significant enhancement can been seen around the transparency 
energy |53l I54j . On the other hand for the steeply falling fluxes the regeneration effect is 
negligible. 



5.2 Neutrino oscillations 

In this section we are going to study how the fluxes of neutrinos at high energies are 
modified by the oscillation mechanism. It is now well established experimentally 55. that 
neutrinos undergo flavour-changing oscillations ■ Due to the high energy of extragalac- 
tic neutrinos, the oscillation length is much larger than the diameter of the Earth so the 
probability that ultrahigh energy neutrinos will oscillate on their way from the surface of 
the Earth to the detector is very small. In fact, the oscillation length in vacuum increases 
linearly with energy 

4ttE . . 

Am; 

If one takes the matter effects into account |57j, then the effective oscillation length in 
matter and the effective mixing angle become 



l-m — lv 



l + 2^cos^ + l f 



„ sin2#„ 

, tan20 m = v , 11 

cos 20 v + l v /L 



where I originates from the matter contribution to the oscillation length and is given by 

V2tt 1.77 x 10 7 
I = T, = m . 12 

G f n e p e 

Here n e is the electron density and p e is the number density in units of Avogadro number/cm 
An important property of l m , is the fact that it saturates at the value I given by (|12[). 
However, at sufficiently high energy l v /l is very large and the effective mixing angle 6 m 
tends to zero. This means that for large energies E > 1 TeV one can safely neglect the 
matter oscillation effects for the neutrinos in their passage through the Earth. 

On the other hand, oscillations can change the initial neutrino flux as compared to the 
primary neutrino flux produced far in the extragalactic source. This is of course due to 
very large distances which neutrinos have to cover on their way from the cosmic source 
to the Earth. It is believed, that the neutrinos emerge as a decay products of the pions 
which are produced in p^ and pp collisions. Since ir^ — > fj, + — > + u e + + 
neutrinos are produced in ratios 1:2:0 for v e : Va : v T . More precisely the v T flux is 
produced [HE] but is smaller than v e and z/„ flux and comes from the decay chain of heavy 
mesons in the pp collisions via the similar process as described in Sec. 2.1 in the case of the 
prompt lepton production in the atmosphere. Due to the very long path to the Earth the 
arriving neutrinos will be completely mixed 59 j which will result in the equal numbers for 
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all species 1:1:1. Let us recall (see for example |59| ) that in the three flavour framework 
the mass and the flavour eigenstates are related by means of the unitary mixing matrix U 











= U 








. u 3 _ 



The probability of the flavour oscillation from the neutrino v a to neutrino v@ after trav- 
elling length L is given by 



P(y a -> u p ;L) = 5, 



" ' _ E U *jUpjU ak U; k (l - exp(-iAE jk L)) . 

J? 4 * 



In the limit when L — > oo one gets 
P(y a -> v p \ L ->• oo) = 5 a p 



l 2 lrr l 2 



(13) 



(14) 



The cosmic neutrino flux F{v a ) after travelling very long distance can be expressed as a 
product of the matrix A and the intrinsic flux F°(u a ) 







' F\v e ) - 




= AA T 




F{v T ) 




F°(u T ) 



where A is defined by (JT 



.4 



Wei? \U e2 \< 



\U, 



If/, 



f*2| 



rl 



r2 



l^3| 2 
r3 



l^3 12 



Since there is an experimental evidence |)U] that ?7 e 3 ~ 0, it means that z/„ and z/ T are 
maximally mixed so that |C/ W -| 2 — | CA- j | 2 < 1 ( j = 1, 2, 3). This symmetry between and 
z^ T results in the complete mixture of the all three neutrino species which would arrive at 
Earth in equal proportion 59_ 1:1:1. 

An interesting possibility of the neutrino decay has been considered recently |61| I62j 
and The two body decays have been considered 



v, 



Vj + X and z/j 



Uj + X , 



(15) 



where X is some very light massless particle, for example a Majoron. It has been shown 
|61| 1621 163j that the decays can lead to ratios which are dramatically different than 1 : 
1 : 1 resulting from oscillations. Thus an important experimental challenge is the flavour 
decomposition of the incoming neutrino flux. The deviation from the uniform mixture 
would definitely point towards new physics like unstable neutrinos or CPT violation. 
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